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1 INTRODUCTION 
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ABSTRACT 

Observations of galaxy luminosity function at high redshifts typically focus on fields of view 
of limited sizes preferentially containing bright sources. These regions possibly are overdense 
and hence biased with respect to the globally averaged regions. Using a semi-analytic model 
based on Choudhury & Ferrara (2006) which is calibrated to match a wide range of observa- 
tions, we study the reionization and thermal history of the universe in overdense regions. The 
main results of our calculation are: (i) Reionization and thermal histories in the biased regions 
are markedly different from the average ones because of enhanced number of sources and 
higher radiative feedback, (ii) The galaxy luminosity function for biased regions is markedly 
different from those corresponding to average ones. In particular, the effect of radiative feed- 
back arising from cosmic reionization is visible at much brighter luminosities, (iii) Because of 
the enhanced radiative feedback within overdense locations, the luminosity function in such 
regions is more sensitive to reionization history than in average regions. The effect of feed- 
back is visible for absolute AB magnitude Mab ^ — 17 at z = 8, almost within the reach 
of present day observations and surely to be probed by the James Webb Space Telescope 
(JWST). This could possibly serve as an additional probe of radiative feedback and hence 
reionization at high redshifts. 

Key words: intergalactic medium cosmology: theory large-scale structure of Universe. 



r e = 0.088 ± 0.015 based on the WMAP seven year data. A 
combination of high redshift luminosity function data with the 
data from these absorption systems and CMB observations favour 
an extende d epoch of reionization tha t begins at 2 « 20 and ends 
at 2 w 6 JChoudhurv & Ferrarall2006h . 



these galaxies, and its evolution, can answer important questions 
about reionization. Indeed, much work has been done on construct- 
ing self-consistent models of structure formation and the evolution 
of ionization and thermal state of the IGM that explain these 
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Studies of the Gunn-Peterson trough ( 


Gunn & Peterson] 096% at 



2^6 have established that the mean n eutral hydrogen fraction is 
higher than 10" 4 (e. g. iFan et alj|2006h and it is most likely that 
the IG M is still highly ionized at these redshifts dGallerani et alj 
l2008al ibh. Furthermore, CMB observations indicate the elec- 
tron scattering optical depth to the last scattering surface to be 
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Nonetheless, interpreting high redshift luminosity functions 
is not straightforward and detailed modelling is required. For in- 
stance, local HII regio ns around these g alaxies can affect luminos- 
ity function evolution JCen et all 2005) and clustering of galaxies 
can enhance this effect JCenll2005h . Another complication is be- 
cause of the fact that these surveys can detect only the brightest 
galaxies at these high redshifts (2 > 6). Such galaxies can form 
only in highly overdense regions and therefore the surveyed volume 
is far from average. An important question in that case is whether 
reion ization proceeds differently in such regions dWvithe & Loebl 
I2OO7I) , 

Galaxy formation is enhanced in overdense regions because 
of a positive bias in abundance of dark matter haloes. The enhance- 
ment in the number of galaxies is proportional to the mass overden- 
sity in the region, with the constant of pro portionality ('bias') re - 
lated to halo masses and collapse redshifts dCoorav & Shethl2002h . 
This increases the number density of sources of ionising radiation 
and aids reionization of the intergalactic medium (IGM) in over- 
dense regions. However, an increase in the IGM density also adds to 
radiative recombinatio n. Furthermore, reionizatio n is accompanied 
by radiative feedback l lThoul & Weinberg! 1 19961) . Radiative feed- 
back heats the IGM and suppresses formation of low mass galax- 
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ies. This increase in radiative recombinations and feedback works 
against the process of reionization and the two effects need not can- 
cel out. Relative significance of these negative and positive contri- 
butions will determine how differently reionization evolves in over- 
dense regions. 

Recently, ItCim et alj d2009l) studied a sample of 4775 -dropout 
candidates identified in five Hubble Advanced Camera for Surveys 
(ACS) fields centred on Sloan Digital Sky Survey (SDSS) QSOs 
at redshifts z ~ 6. They compared results with those from equally 
deep Great Observatory Origins Deep Survey (GOODS) observa- 
tions of the same fields in order to find an enhancement or suppres- 
sion in source counts in ACS fields. An enhancement would imply 
that bias wins over negative feedback in these overdense regions. 
They found the ACS populations to be overdense in two fields, un- 
derdense in two field, and equally dense as the GOODS popula- 
tions in one field. Somewhat surprisingly, they did not find a clear 
correlation between density of 2775 dropouts and the region's over- 
density. 

In this paper, we use semi-analytic models to study reioniza- 
tion within overdense regions. The main aim of this work is to 
quantify the effects of enhancement in the number of sources and 
radiative feedback within such regions and explore the possibility 
whether the galaxy luminosity function in overdense regions can 
be used as potential probe of feedback and reionization history. It is 
known that if ionization feedback is the main contributor to the sup- 
pression of star formation in low mass haloes then one can distin- 
guish between early and late reionization histories by constraining 
the epoch at which feedback-related low-luminosity flattening oc- 
curs in the galactic luminosity function. The effect of reionization 
feedback on the high redshift galaxy luminos ity function was firs t 
demonstrated using semi-analytic models by ISamuiet "all d2007l) . 
We apply their method to study the luminosity function in over- 
dense regions. 

We describe our model for reionization and give details about 
various parameters and their calibration in §2. We explain our 
method of identifying biased regions and outline modification to 
our reionization model in such regions in §3. Details of our lumi- 
nosity function calculation appear in §4. We present and discuss our 
results in §5 and §6. Throughout the paper, we use th e best-fit cos- 
mological parameters from the 7-year WMAP data dLarson et al.l 
l20ld) . i.e., a flat universe with Q, m = 0.26, tt A = 0.73, Q K = 
0.044 and Q b = 0.045, and h = 0.713. The parameters defining 
the linear dark matter power spectrum are as = 0.80, n s = 0.96, 
dn s /dlnfc = -0.034. 



2 DESCRIPTION OF THE SEMI-ANALYTIC MODEL 

In this section, we first summarise the basic features of the semi- 
analytic model used for studying the globally averaged reionization 
history. We then describe in detail the modifications made to this 
model in order to study reionization in biased regions. 



2.1 Globally averaged reionization 

Our model for reionization an d thermal history of the avera ge IGM 
is essentially that developed in lChoudhurv & Ferraral2005l (CF05). 
The main features of this model are as follows. 

The model accounts for IGM inhomogeneities by adopting a 
lognormal distribution with the evolution of volume filling factor of 
ionized hydrogen (Hll) regions Qhii(z) being calculated according 



to the method outlined in Miralda-Escude et al. I feOOOh : reioniza- 
tion is said to be complete once all the low-density regions (say, 
with overdensities A < A cr i t ~ 60) are ionised. We follow the 
ionization and thermal histories of neutral and HII regions simul- 
taneously and self-consistently, treating the IGM as a multi-phase 
medium. In this work, we do not consider the reionization of singly 
ionised helium as it occurs much later (z ~ 3) than redshifts of our 
interest. 

The number of ionising photons depends on the assumptions 
made regarding the sources. In this work, we have assumed that 
reionization of hydrogen is driven by stellar sources. The rate of 
ionising photons injected into the IGM per unit time per unit vol- 
ume at redshift z is denoted by hph(z) and is essentially deter- 
mined by the star formation rate (SFR) density p,(z). The first 
step in this calculation is to evaluate the comoving number density 
N(M, z, z c )dMdZc at redshift z of collapsed halos having mass in 
the range M a nd M + dM and redshift of collapse in the range z c 
and z c + dz c dSasakilll994h : 



N(M, z, z c )dMdz c = N(M, z c )v(M, z c ) 



D(z c ) 
D(z c ) 



x psurv(z, z c )-j^-dz c dM, (1) 

dZc 

where N(M, z c )dM is the comoving number density of collapsed 
halos with mass between M and M + dM, also known as the 
Press-Schechter (PS) mass function dPress & Schechte3ll974h . and 
p S urv(z, z c ) is the probability of a halo collapsed at redshift z c sur- 
viving without merger till redshift z. This survival probability is 
simply given by 



rv Zc) 



D(z) 



(2) 



where D(z) is growth function of matter perturbations. Further- 
more, v(M, z c ) is given by S c /[D(z c )a(M)], where a(M) is the 
rms value of density fluctuations at the comoving scale correspond- 
ing to mass M and 8 C is the critical overdensity for collapse of the 
halo. Next, we assume that the SFR of a halo of mass M that has 
collapsed at an earlier redshift z c peaks around a dynamical time- 
scale of the halo and has the form 



x exp 



t(z)-t(z c ) 



^dyn(-2c) 



(3) 



where /* denotes the fraction of the total baryonic mass of the halo 
that gets converted into stars. The global SFR density at redshift z 
is then 



P*(z) = 



dz c 



dMM*(M,z,Zc)N(M,z,z c ), (4) 



Af min ( 2c ) 



where the lower limit of the mass integral, A/ m i n (z c ), prohibits 
low-mass halos from forming stars; its value is decided by differ- 
ent feedback processes. In this work, we exclusively consider ra- 
diative feedback. For neutral regions, we assume that this quantity 
is determined by atomic cooling of gas within haloes (we neglect 
cooling via molecular hydrogen). Within ionised regions, photo- 
heating of the gas can result in a further suppression of star for- 
mation in low-mass haloes. We compute such (radiative) feedback 
self-consistently from the evolution of the thermal properties of the 
IGM, as discussed in Section [23l 

We can then write the rate of emission of ionising photons per 




Figure 1. Behaviour of various quantities in our fiducial model, in the average and overdense reg ions are shown by dashed and solid lines respectively. The top 
left panel shows the photoionisation rate, with data points taken from IBolton & Haehnelt 2007). The top right panel shows the mass-averaged temperature for 
ionised regions, which essentially determine the radiative feedback. The bottom left panel is for the volume filling factor of ionised regions. The bottom right 
panel shows the cosmic star formation rate. Note that the overdense region that we consider here collapses at z = 6.8. The vertical dotted line in the top left 
panel highlights this. We cannot evolve our reionization model for smaller redshifts. 



unit time per unit volume per unit frequency range, n„(z), as 

h v (z) = 7Y 7 (y)/escP*(^) ; (5) 

where N-,(u) is the total number of ionising photons emitted per 
unit frequency range per unit stellar mass and / csc is the escape 
fraction of photons from the halo. The quantity N 1 {p) can be 
calculated using population synthesis, given t he initial mass func - 
tion and spectra of stars of different masses JSamui et al.ll2007h . 
In this paper we have used t he population synthesis co de Star- 
burst ^ dLeitherer et alll 19991 : rvizquez & Leithererll2005l) to cal- 
culate N 1 {v) by evolving a stellar population of total mass 10 6 
M with a 0.1 - 100.0 M Salpeter IMF and metallicity 0.001 
(0.05 times the solar metallicity, Zq — 0.02). The total rate of 
emission of ionising photons per unit time per unit volume is ob- 
tained simply integrating by Equation (O over suitable frequency 
range. 

Given the above model, we obtain best-fit parameters by com- 
paring with the redshi ft evolution of photoionisa tion rate obtained 
from the Lya forest ( IBolton & H aehnelt 2007) and the electron 
scattering optical depth I Larson et al. 2010h . We should mention 
here that any model containing only a single population of atomic- 
cooled stellar sources with non-evolving .f t/esc cannot match both 
the Lya forest and WM AP constraints dChoudhurv et alj 120081 ; 
IBolton & Haehneltl2007h . In this work, we choose the model which 
satisfies the Lya constraints but underpredicts r e . In order to match 
both the constraints, one has to invoke either molecular cooling in 



minihaloes and/or metal-free stars and/or other unknown sources 
of reionization. This model is described by the parameter values 
/» = 0.2 and / esc = 0.135, and gives r e of 0.072. Figure Q] shows 
evolution of the filling factor of ionised regions, global star for- 
mation rate density, mass-weighted average temperature in ionised 
regions and average hydrogen photoionisation rate in this model 
(dashed curves in all the panels). The filling factor of ionized re- 
gions is seen to rise monotonically from z ss 15 and takes values 
close to unity at redshifts z ~ 6. Temperature of ionized regions 
also rises rapidly during reionization and flattens out to a few times 
10 4 K at redshift z < 4 (not shown here). Lastly, the photoioni- 
sation rate also increases during reionization as the star formation 
rate builds up. However, the photoionisation rate increases rapidly 
with a sudden jump at z ~ 6 when the ionized regions overlap (fill- 
ing factor becomes close to unity). This is because a given region 
in space starts receiving ionizing photos from multiple sources and 
as a result, the ionizing flux suddenly increases. This is our fidu- 
cial model, which satisfies observational constraints from Lya for- 
est, observations of star formation rate history, number density of 
Lyman-limit systems at high redshift and of the IGM temperature. 
In this model, reionization starts at z ~ 15 and is 90% complete by 
z ~ 7. Evolution of xhii is consistent with constraints from Lya 
emitters and the GP optical depths. 

Having set up the reionization model, we then calculate the 
predicted luminosity function of galaxies in this model. Luminosity 
functions of objects are usually preferred for comparing theory with 
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observations because of its directly observable n ature. In thi s work , 
we closely follow the approach presented by ISamui et alj J2007h 
to calculate the luminosity function. We obtain luminosity per unit 
mass, lisao(t), at 1500 A as a function of time from population 
synthesis for an instantaneous burst. In our model, star formation 
does not happen in a burst, but is a continuous process spread out 
over a dynamical time-scale. Therefore, in order to determine the 
luminosity of a halo, L1500 (t) with this kind of star formation, we 
convolve £1500 (t) with the halo's star formation rate using 



(6) 



L 1500 (M,T)= / dTM„(M,T -t,z c )Ii5oo(t), 



where T is the age of the halo, which has mass M and which col- 
lapsed at redshift z c . This luminosity can be converted to absolute 
AB magnitude using 



Mai 



-2.5 log 10 (L„ ) + 51.60, 



(7) 



where the luminosity is in units of erg s _1 Hz -1 dOke & GunrJ 
1983). One can compute the luminosity evolution for any halo that 
collapses at redshift z c and undergoes star formation according to 
Equation l[3}. The luminosity function at redshift 2, $(JWab, 2), is 
now given by 



<f>(MAB,z)dM A B 

dz c N(M, 2, z c 



dM dLi 5 



dh x 



dM Aj 



-dM AB , 



(8) 



where N(M, z, z c ) is the number density at redshift 2 of halos of 
mass M collapsed at redshift z c . We will use Equation {8} to study 
effect of overdensity on the luminosity function and to compare 
the luminosity function in our model with observations in the next 
section. 

Figure[2]shows the globally averaged luminosity function cal- 
culated using our model fo r different redshifts in compar ison with 
observations presented bv lBouwens & Illingworth (2006). We find 
that our model reproduces the observed luminosity functions at 
high redshifts reasonably well. In particular, the match at 2 = 6 is 
remarkably good while the model predicts less number of galaxies 
than what is observed at 2 = 7 and 8. This could indicate that the 
star-forming efficiency /„ increases with 2, and/or the time-scale of 
star formation is lower than td yn at higher redshifts. The match of 
the model with the data can be improved by tuning these parame- 
ters suitably, however we prefer not to introduce additional freedom 
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Figure 3. Effect of overdensity on the luminosity function via feedback at 
z=8.0. This is for <5 = 8.8 and R L = 1.482 Mpc. 



in constraining the parameters; rather our focus is to estimate the 
effect of reionization and feedback on the luminosity function. 

In our calculation of luminosities, we do not make any correc- 
tion for dust. This is partly beca use of indications from observed 
very blue UV-continuum slopes iBouwens et alj|2010l : lOesch et"al] 
l20ld : lFinkelstein et alJl20ld:lBunker et alj|20irj|) that dust extinc- 
tion in 2 > 7 is small. As discussed in the next section, we exclu- 
sively work with luminosity functions at these redshifts. Also, the 
effect of dust is degenerate with /« to some extent. Therefore, the 
exclusion of dust extinction does not affect the general results of 
our calculation. 

2.2 Biased regions 

We have mentioned that galaxy luminosity functions provide valu- 
able information regarding reionization. However, observations are 
carried out over relatively small fields of view. The bright sources 
in these fields are typically hosted by high mass haloes. Hence, 
it is most likely that these fields are biased tracers of luminosity 
function and reionization. In this section, we extend our model to 
study reionization within such biased regions and quantify the de- 
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parture of various quantities from their globally averaged trends. 
Reionization in biased regions has been discussed in the literature. 
Wvith e & Loebl d2007h studied the correlation between high red- 
shift galaxy distribution and the neutral Hydrogen 21 cm emis- 
sio n by considering re ionization in the vicinity of these galax- 
ies, considered the ioni zation background near 
high redshift quasars. lGeil& Wvithe] 12009) studied effect of reion- 
ization around high re dshift quasars on the power s pectrum of 21 
cm emission (see also IPritchard & Furlanettd 120071) . The general 
conclusion of these studies is that overdense regions are ionised 
earlier. In this work, we will consider the behaviour of luminosity 
functions in such regions. 

Overdense regions are characterised by their comoving La- 
grangian size R, and their linearly extrapolated overdensity S. At a 
given redshift, we can take a scale R corresponding to an observed 
field of view (e.g., WFC3/IR field in HST) and then determine 5 
by identifying the presence of a massive object, for example, a 
quasar or a bright gal axy. We follow a prescription discussed by 
iMunoz & LoerT l l20()8h . 

Note that if a galaxy with luminosity L1500 is observed at red- 
shift z g then we can assign a certain mass to the dark matter halo 
containing the galaxy, say M. The halo mass M has to be obtained 
from galaxy luminosity L1500 by using some pres cription or by fit- 
ting th e galaxy's spectral energy distribution (SED; IVale & Ostrikeil 
2004). In this work, however, since uncertainty in the value of over- 
density 8 is expected to be larger than any uncertainty in the mass of 
the galaxy's host halo, we choose to calculate this mass in an alter- 
nate, simpler manner. Note that we can invert equations ((6} and {3} 
to obtain M given an observed value of L1500 if we have an esti- 
mate of the redshift at which the galaxy formed. In fact, we obtain 
the halo mass M such that the halo luminosity after a dynamical 
time from halo formation time, calculated according to our model, 
is equal to Lisoo • In other words, we break the degeneracy between 
halo mass and formation redshift by assuming that the galaxy's age 
is equal to the dynamical time of the halo. 

From this analysis, we conclude that a collapsed object with 
mass M exists at redshift z g . Now suppose our field of observa- 
tion corresponds to some linear scale R at this redshift. Then the 
linearly extrapolated overdensity at this scale can be obta ined us- 
ing the excursion set prescription flMunoz & Loebl 120081) . Recall 
that the probability distribution of the extrapolated Gaussian den- 
sity field smoothed over scale R is also Gaussian 



Q{8 R ,a(R))dS R 



■. exp 



2<r 2 (7?) 



d8 f 



(9) 



The conditional probability distribution of overdensity 8\ on a scale 
specified by the variance a\, given a value of overdensity 82 on a 
larger scale specified by the variance <x§ < g\ is given by 



Q(8\,a\\8i-,o\) = Q(Si - S 2 , 



2 2n 
<? 1 - 0a) 



(10) 



Conversely, when the value of overdensity Si on a smaller scale 
specified by variance a\ is given, the conditional probability distri- 
bution of 82 can be obtained using Bayes theorem as 



Q(<52,02|<5l,£Ti) OC Q (<5l,(7l 1 82,02 )Q ($2,02 )dS2. 



(11) 



If we now set the smaller scale to be that of the observed collapsed 
halo, and the overdensity at that scale to be the critical overden- 
sity for spherical collapse, Equation il lb will give the resulting 
overdensity at any larger scale due to the presence of this massive 
galaxy. In other words, we set 81 = 8 c (z g ) and a\ = <j 2 (M) in 
Equation dl It . 

The larger scale corresponds to the field of observation. In 



order to calculate that, we first note that the excursion set prin- 
ciple functions entirely in Lagrangian coordinates. As a region 
evolves towards eventual collapse its Lagrangian size stays un- 
changed while its Eulerian size changes. For a spherical region the 
Eulerian evolution will follow the solution of the spherical collapse 
model. However, since the Eulerian and Lagrangian sizes of the re- 
gion coincide at the initial instant, the spherical collapse solution is 
also a relationship between these two sizes. Thus we have 



Re = 



3 1- cos 6D{z = 0) 



10 



D{z) 



Rl, 



where 8 is a parameter given by 

1 _ 3x6 2 / 3 (0-sin0) 2 / 3 



1 + z 



20 



(12) 



(13) 



In our case, the Eulerian size of the region of interest is just the 
comoving distance corresponding to the angular field of view, 
which is just the angular diameter distance at the relevant red- 
shift multiplied by the angular field of view. The WFC3/IR field 
is 136 x 123 . For the best fit ACDM cosmology the diago- 
nal size of this field corresponds to a comoving Eulerian distance 
Re = 1-365 Mpc at z = 8. In a WFC3/IR field centred on 
the object UDFy-42886345 at redshift 8.0 and apparent magnitude 
Hieo,AB = 28.0 we obtain a halo mass M = 2.52 x 10 11 M 
and luminosity L1500 = 2.21 x 10 29 erg s _1 Hz -1 . 

Notice, however, that since 8l is unknown, Equation j 1 2t im- 
plies that the relation between the Eulerian size Re and Lagrangian 
size Rl is not one-to-one. Thus, for the probability distribution of 
linearly extrapolated overdensity 8 given the halo mass M, we can 
only write 



dP(8\M) 
d8 



<xQ[S,R L (S, Re)\S c (z),R(M)] 



(14) 



where the constant of proportionality is calculated by using the nor- 
malization condition J[dP(S\M)/dS]dS — 1. 

In our calculations, we work with the value of 8 for which 
dP(S\M)/dS is maximum. For the WFC3/IR field at z = 8, this 
turns out to be 8 = 8.86 (linearly extrapolated to z = 0), which 
results in a Lagrangian size Rl ~ 1.482 Mpc for the region of 
interest. Notice that since 8 > S c the region must have collapsed at 
some redshift z < 6.5. 

In order to incorporate this overdensity into our reionization 
model, note first that the number density of collapsed objects in 
such overdense regions is enhanced with respect to that in a region 
with average density. Thi s enhancement ca n be calculated using the 
excursion set formalism dBond et al ] |l99ll) . It is then shown in the 
Appendix that the comoving number density N(M, z, z c )dMdz c 
at redshift z of collapsed halos having mass in the range M and 
M + dM and redshift of collapse in the range z c and z c + dz c is 
given in this case by 



N(M, z, z c )dMdz c = N(M, z c ) 



v 2 8 c 



D(zc) 
S c /D(z c )-Sj D 2 (z c ) 
dt 



x Psurv(«, z c )-—dz c dM, 
dz c 



(15) 



where N(M, z c ) is the PS mass function. The scale R enters via 
the definition of v(M, z c ), which is now given by 



v{M, Zc) 



8 C /D{z c ) - 8 



(16) 
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The survival probability p surv ( 
S c /D(z) - 8 



is given by 



Psurv ('2- 2 C ) 



S c /D(z c ) — S 



(17) 



Another change when our reionization model is applied to 
overdense regions is that we now normalize the probability distri- 
bution of inhomogeneities in the IGM such that the average density 
in the region is A = 8 + 1. 



2.3 Radiative feedback 

As we argue in the next section, the luminosity function of galaxies 
in an overdense region could carry an enhanced signature of feed- 
back. We therefore highlight our feedback model in this subsection. 

Radiation from stars in the first galaxies is expected to ionize 
and heat the surrounding medium. This increases the mass scale 
above which baryons can collapse in haloes. Also, as a result, the 
minimum mass of haloes that are able to cool is much higher in ion- 
ized regions than in the neutral ones. In our calculations, feedback 
appears through the quantity M m i n (z) in Equation|4] The tempera- 
ture evolution of both regions is calculated self-consistently. In the 
ionized regions, we fix the cut-off mass to that corresponding to a 
virial temperature of 10 4 K or the local Jeans mass, whichever is 
higher. In the neutral regions, since the Jeans mass is always low, 
the cut-off mass always corresponds to the virial temperature of 
10 4 K. The minimum mass corresponds to the circular velocity of 



2/Cbo 



/j,m p 



(18) 



where /i is the mean molecular weight. For a temperature of ~ 
10 4 K, the minimum circular velocity is ~ 25 km s . Note that 
this v alue is comparable to values obtained in simulations ( Gnedin 
200G) but is s omewhat higher tha n that taken in the semi-analytic 
prescription of lSamui et alJ d2007h . 

We find that M m i n (z) increases with time taking values of 
fa 10 7 M at z w 10 and » 10 s M at z fa 7. In overdense 
regions the minimum mass is enhanced to about 10 10 Mq. Figure 
f5]shows the evolution of the minimum mass. 



3 RESULTS 

The results for reionization and thermal histories within overdense 
regions are presented in this section. 



3.1 Effect of overdensity on reionization history 

We first consider the effect of overdensity on reionization history 
for our fiducial model. As is well known, reionization proceeds 
differently in overdense regions. The solid lines in Figure [T] show 
the evolution of the photoionisation rate, temperature in ionised re- 
gions, star formation rate density and the volume filling factor of 
ionised regions in an overdense region with size Rl — 1.482 Mpc 
and linearly extrapolated overdensity 8 = 8.86. This corresponds 
to the HUDF WFC3 /IR field centred at the brightest source in 
iBouwensetaflfcOld) . (See Section |2~2l ) Clearly while the average 
region is completely ionised at z w 6, the biased region is ionised 
much earlier, at z w 7.5. This result agrees with Wvit he & Loebl 
although note that unlike that work, here we calculate the 
clumping factor from a physical model for inhomogeneities. The 
reason for early reionization in overdense regions is the enhanced 
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Figure 5. Evolution of the minimum mass M m i n (z) of haloes that can 
host galaxies in ionized regions in the average case (dashed line) and the 
overdense case (solid line). The dot-dashed line shows the minimum mass 
in neutral regions, which is same for the average and overdense cases. 



number of sources, which is clear from the plots of photoionisa- 
tion rate and the star-formation rate, both of which are ~ 5 times 
higher than the corresponding globally averaged values. However, 
these overdense regions have more recombinations, which results 
in enhanced temperatures as is clear from the top right panel. This 
results in enhanced negative radiative feedback which will suppress 
star formation in low mass galaxies and hence affect the shape 
of the luminosity function. In fact, for the average case, haloes 
in ionised regions with masses below 10 8 Mq cannot form stars, 
whereas this cutoff mass rises to close to 10 10 Mq in the overdense 
case. Clearly feedback is enhanced in overdense regions. 

3.2 Effect of overdensity on luminosity function 

We now discuss the effect of overdensity on luminosity func- 
tion. Clearly, overdense regions have enhanced number of sources, 
hence it is natural that the amplitude of the luminosity function for 
such regions should be higher than the globally averaged values. 
However, the overdense regions have enhanced radiative feedback 
too, and hence we expect a decrease in the number of sources, par- 
ticularly towards the fainter end. 

Figure [3] shows the effect of overdensity on the luminosity 
function at z — 8 for our fiducial model. The ionised volume filling 
factor within the overdense region is Q fa 1.0 for the overdense re- 
gion under consideration at this redshift. The average region lumi- 
nosity function (dashed line) is clearly very different from the lumi- 
nosity function in the overdense region (solid line) at that redshift. 
Firstly, we can clearly see an enhancement in the source counts for 
brighter galaxies, which is as expected. In addition, there is a clear 
sign of a flattening for magnitudes Mab > —17, which is a signa- 
ture of radiative feedback. In comparison, the effect of feedback for 
average regions occurs at much fainter magnitudes Mab ~ —12. 
Note that there is no complete suppression of star formation for 
halo masses lower than the feedback threshold, rather the luminos- 
ity function for magnitudes below the knee continues to grow in the 
flattened region. This is simply due to the continued star formation 
in haloes with mass less than the cutoff mass at z = 8.0, but which 
collapsed at higher redshifts when the feedback threshold mass was 
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Figure 4. Effect of reionization history on luminosity function at z=8.0. Solid, dashed and dot-dashed lines have r = 0.073, 0.058 and 0.088 respectively. Left 
panel shows the average case. Right panel shows the overdense case. 



lower. Thus, for instance, if star formation is allowed to happen in a 
halo for only for a fraction of the dynamical time [see equation l[3}], 
the luminosity function will rise less steeply at the fainter end. For 
small enough star formation time scale, the luminosity function will 
show an abrupt cutoff. Of course, an abrupt cutoff is always seen at 
low enough luminosities, which are not shown in the figure here. 

It is important to understand here that the data points in Figure 
[2]do not represent luminosity function of the overdense region. In- 
stead, those data points represent the globally averaged luminosity 
function derived using a maximum likelihood procedure from the 
observed luminosity distribution of sources. In this procedure, a 
likelihood function is defined, which describes the step-wise shape 
of the luminosity function that is most likely given the observed lu- 
minosity distribution in the search fields. D etails of this procedure 
are described, for example, in Section 5.1 of lBouwens et alj J201C|) 
and references therein. 



3.3 Luminosity function as a probe of reionization 

Given the fact that the effect of radiative feedback shows up at 
brighter magnitudes for overdense regions, it is possible to use this 
feature for studying feedback using near- future observations. For 
this purpose, we consider two additional models (other than the 
fiducial one) of reionization. These models have parameter values 
(/», /esc) = (0.06, 0.3) and (0.2, 0.07) and we obtain r e = 0.088 
and 0.058 respectively for these models. We fix /, and only change 
the value of / csc to ensure that any effect on the luminosity function 
is purely due to feedback. These two models predict photoionisa- 
tion rates greater and lesse r respectively than what are presented by 
iBolton & Haehnelj fcOOTh . 

The right panel of Figure [4] shows the luminosity function at 
z = 8 within the overdense region for three different reionization 
histories, which can be compared with the corresponding luminos- 
ity function in average region (shown in the left panel). In both 
cases a distinct "knee" is seen in the luminosity function as a signa- 
ture of feedback. The luminosity function flattens at this luminosity, 
and is suppressed to very low values at much lower luminosities. As 
described in the previous section, this signature of feedback appears 



at brighter magnitudes for the overdense region. This is expected, 
because the cutoff mass depends directly on the temperature, which 
is enhanced in the overdense region. We also note that in the case 
of the first model the flattening occurs for Mab > —19 whereas 
for the second model at a fainter luminosity of Mab — —16. This 
is due to the fact that the photoionisation feedback is enhanced in 
the first model due to enhanced flux. 

The evolution of the filling factor affects this result through 
the average temperature which sets the cutoff mass. Thus, early 
and late reionization models are distinguished by the difference in 
the nature of flattening in both cases. This also affects the evolution 
of the luminosity function. 

We find that the reionization history has a strong effect on the 
luminosity function at the faint end. It is known that the bright end 
of the luminosity function is affected primarily by the star forma- 
tion mode of a halo, and the overall bias, whereas its faint end is 
affected by the reionization history. 

However, we also find, from Figure[4] that the effect of reion- 
ization history is much stronger in the case of overdense regions. 
This is because of the enhanced photoionisation feedback, which 
is more sensitive to changes in reionization history. This order 
of magnitude change in the overdense region luminosity function 
should be visible to the James Webb Space Telescope, which can 
observe up to ttlab ~ 31.5 (Mab ~ —16.0 atredshifts of interest; 
IWindhorst et"al]2006h . 



4 DISCUSSION AND SUMMARY 

We have used a se mi-ana lytic model, based on 

Choudhur v~& Ferraral j2005l [2006) to study reionization and 
thermal history of an overdense region. Studying such regions 
is important because observations of galaxy luminosity function 
at high redshifts typically focus fields of view of limited sizes 
preferentially containing bright sources; these regions possibly are 
overdense and hence biased with respect to the globally averaged 
regions. In particular, we study the effect of radiative feedback 
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arising from reionization on the shape of galaxy luminosity 
function. 

In summary, we find that 

(i) Reionization proceeds differently in overdense regions. 
Overdense regions are ionised earlier because of enhanced number 
of sources and star formation. In addition, these regions have higher 
temperatures because of enhanced recombinations and hence effect 
of radiative feedback is enhanced too. 

(ii) In particular, the shape of the galaxy luminosity function for 
biased regions is very different from that for average regions. There 
is a significant enhancement in the number of high-mass galaxies 
because of bias, while there is a reduction in low-mass galaxies 
resulting from enhanced radiative feedback. 

(iii) Luminosity function in overdense regions is more sensitive 
to reionization history compared to average regions. The effect of 
radiative feedback shows up at absolute AB magnitudes Mab > 
— 17 in these regions, while it occurs at much fainter magnitudes 
Mab ~ — 12 for average regions. This order of magnitude change 
in the overdense region luminosity function should be visible to the 
James Webb Space Telescope in future. 

Finally, we critically examine some of the simplifying as- 
sumptions made in this work and how they are likely to affect our 
conclusions. Firstly, we have seen that the presence of a high mass 
galaxy within a region if size R does not uniquely specify the value 
of the overdensity 8. Rather we obtain a probability density (which 
is Gaussian in shape) and work with the value where this proba- 
bility is maximum. In reality, however, the actual value of 5 could 
be different and this may possibly affect the predicted luminosity 
function. Note that the luminosity function at the brighter end is al- 
most independent of the details of reionization history, and this, in 
principle, can be used for constraining the value of 8. The effect of 
feedback can then be studied using the faint end of the luminosity 
function. 

The radiative feedback pr escription used in this paper is based 
on a Jeans mass calculation ( Choud hury & Ferrar al l2005h . How- 
ever, al t ernate prescriptions for feedback exist in literature, e.g., 
iGnedinl fioooh and hence the shape of the luminosity function at 
faint ends as predicted by our model may not be robust. Interest- 
ingly, the presence of a "knee" in the luminosity function can be 
used to estimate the value of the halo mass below which star forma- 
tion can be suppressed (which in turn can indicate the temperature) 
while the shape of the function below this knee should indicate 
the nature of feedback. This study can also be complemented with 
proposed for st udying feedback using other obser vations, e.g., 21 
cm ob servation {Schneider et al. 2008) and CMBR dBurigana et al.l 
l2008h . 

Finally, we have neglected the presence of other sources of 
reionization, e.g., metal-free stars, minihaloes, and so on. It is ex- 
pected that these sources would be too faint to affect the luminosity 
function in the ranges we are considering. However, these sources 
may affect the thermal history of the medium, e.g, the metal-free 
stars would produce higher temperatures because of harder spec- 
tra. In such cases, it is most likely that feedback would occur at 
magnitude brighter than what we have indicated and hence would 
possibly be easier to detect. 
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APPENDIX: FORMATION RATE AND SURVIVAL 
PROBABILITY OF HALOES IN OVERDENSE REGIONS 

As expressed in Equation {5}, the number density of ionizing pho- 
tons produced per unit time is related to the SFR density, which 
in turn depends on the SFR in each halo, given by Equation (01, 
and the number density of haloes of a certain age, given by Equa- 
tion Q} for average regions, and by Equation dl5t for overdense 
regions. We derive Equation J15l > in this appendix. 

We denote the number density at redshift z of haloes formed 
between redshifts z c and z c + dz c , with mass between M and 
M +dM, by N(M, z, z c )dMdz c . This quantity is related to (1) the 
formation rate at redshift z c of haloes with mass between M and 
M + dM, denoted by iV form (M, z c )dM, and (2) the probability 
of their survival at redshift z, denoted by p SUIV (z, z c ). We calcu - 
late these two quantities using a technique given bv lSasakilfl994» . 
applied to an overdense region with overdensity S and size R. 

Recall that in extended Press-Schechter theory dBond et all 

Il99lh . the mass function of dark matter haloes is defined as the 
comoving number density of haloes with mass between M and 
M + dM. At redshift z, this quantity is given by 

N(M, z)d M = ^ex P (z£)^ d M, (19) 

where p m is the average matter density, and, as before, v{M, z) = 
S c /[D(z)a(M)]. The critical overdensity of collapse of a halo is 
denoted by 8 C , D(z) is the growth function of density perturba- 
tions, and <r(M) is the rms value of density perturbations at the 
comoving scale corresponding to mass M. In a region with over- 
density S and linear size R, the mass function is enhanced. This 
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enhancement can be calculated using the excursion set formalism 
jBond et al.lfl99lh . The resulting mass function is again given by 
Equation d!9t . except that now the quantity v{M, z) is defined as 



Sc/D(z 



(20) 



V^(M) - o\ 

where or is the rms valu e of density p erturbations at comoving 
scale R. Closely following lSas we can write 

N(M, z) = N lom (M, z) - N dcst (M, z), (21) 

where Nd BB t(M, z)dM is the destruction rate at redshift z of haloes 
of mass between M and dM. (The halo formation rate is defined as 
the number density of haloes formed per unit time from mergers of 
lower mass haloes. Similarly the halo destruction rate is defined as 
the number density of haloes destroyed per unit time due to mergers 
with other haloes.) Here, an overdot denotes the time derivative. We 
can write the destruction rate as 



N dest (M, Z ) = 



N(M,z)Q(M,M',z)dM', 



<j)(M,z)N(M,z) 



and the formation rate as 

Nf om (M,z)= / N{M',z)Q(M',M,z)dM', 



(22) 
(23) 

(24) 



where Q(M, M' , z) is the probability that a halo of mass M 
merges with another halo to result in a halo of mass M' per unit 
time, and Q(M' , M, z) that an halo of mass M forming at red- 
shift z has a progenitor of mass M' . The threshold mass M min is 
introduced at this stage to avoid divergence. This gives 



7V form (M, z) = N(M, z) + <j>{M, z)N(M, z). 



(25) 



We now assume that <j> has no characteristic mass scale so that 
<(>{M, z) = M a <f>(z). This gives 

7/ * -N(M, z) + N lmm (M,z) 

^ = N(M,z)M" • (26) 

But since the left hand side of Equation l |26l l is a function of time 
alone (through the redshift), the right hand side of this equation 
also has to be independent of mass. In particular, we can then set 
M — Af m i n in this equation, giving us 



-jV(M min , z) 



(27) 



^ ' N(M min ,z)M« in - 
Now, in the case of the overdense region that we are considering 



here, we have 



D(z) 



D 2 (z) 8 C /D(z)~8 



[v 2 {m,z)-l], (28) 



[v 2 (M min ,z)-l]M- 



(29) 



N(M,z) = N(M,z] 

which gives 

7_ D 8 C 
D 2 S c /D(z)-5 l 

Since our choice of threshold mass M m i n is arbitrary, we now need 
to take the limit M m m — > 0. However, since v — > in this limit, 
<f> becomes indeterminate, except when a — 0. This implies that 
we must set a — for consistency. This gives (j>(M, z) = 4>{z). 
Substituting the resultant expression in Equation l |25t , we get 



N tom {M,z) = N(M,z) 



D 



5c 



-.u{M,z). 



(30) 



D 2 8 C /D(z) - S 
This the required formation rate of haloes in an overdense region 



Furthermore, from our definitions of probabilities in Equa- 
tions J23b and l !24t . we can write the probability that a halo that 
has formed at redshift z c continues to exist at redshift z as 



,{z, z c ) = exp 



t(z) 



<j>(t')dt' 



't(Zc) 

which in our case results in 

S c /D(z) - S 



(3D 



(32) 



S c /D(z c )-S- 

From Equations J30t and J32b . we can now write the the comov- 
ing number density N(M, z, z c )dMdz c at redshift z of collapsed 
halos having mass in the range M and M + dM and redshift of 
collapse in the range z c . and z c + dz c as 



N(M, z, z c )dMdz c = N(M, Zc) 
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This is our Equation (15}. 

It is worth pointing out that Equations d30t and J3 2. b reduce to 
the average forms for halo formation rate and survival probability 
in the limit 6 — > and R — > oo. 



